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Abstract 

We propose a simple geometrical approach for finding the Lewenstein-Sanpera decomposition 
of Bell decomposable states of 2 <g> 2 quantum systems. We show that in these systems, the 
weight of the pure entangled part in the decomposition is equal to the concurrence of the state. 
It is also shown that the optimized separable part of L-S decomposition minimizes the von 
Neumann relative entropy. We also obtain the decomposition for a class of mixed states by 
using some LQCC actions. It is also shown that for these states the average concurrence of L-S 
decomposition is equal to their concurrence. 

Keywords: Quantum entanglement, Bell decomposable states, Lewenstein-Sanpera 
decomposition, Concurrence 
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1 Introduction 

Entanglement is one of the most striking features of quantum mechanics [Q, |2| . The non local char- 
acter of an entangled system is usually manifested in quantum correlations between non interacting 
subsystems provided that they had only interaction in the past. A bipartite mixed state is said to 
be separable (non entangled) if it can be written as a convex combination of pure states 



where 



and 



iftf/ are pure states of subsystems A and B, respectively. In the case of pure states 
it is easy to check whether a given state is, or is not entangled. Entangled pure states do always 
violate Bell inequalities 0. For mixed states, however, the statistical properties of the mixture 
can hide the quantum correlations embodied in the system, making thus the distinction between 
separable and entangled states enormously difficult. 

In the pioneering parer H, a very interesting description of entanglement was achieved by 
defining the best separable approximation (BSA) of a mixed state. In the case of 2-qubit system, 
it consists of a decomposition of the state into a linear combination of mixed separable part and a 
pure entangled one. In this way, the whole non-separability properties are concentrated in the pure 
part. It also provides a natural measure of entanglement given by the entanglement of the pure 
part (well defined for pure states) multiplied by the weight of the pure part in the composition. 

In the Ref. the numerical method for finding the BSA has been reported. Also in 2 (g) 2 
systems some analytical results for special states were found in Q. An attempt to generalize the 
results of Ref. [Q] is made in || . 

In an algebraic approach to find BSA of a 2-qubit state is attempted. They have also showed 
that the weight of the entangled part in the decomposition is equal to the concurrence of the state. 
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In this paper we consider Bell decomposable (BD) states. We provide a simple geometrical ap- 
proach and give an analytical expression for L-S decomposition, where our results are in agreement 
with those reported in [Q, || . Our method to find L-S decomposition is geometrically intuitive. We 
also see that the weight of the entangled part in the decomposition is equal to the concurrence 
of the state. It is also shown that separable state optimizing L-S decomposition, minimizes the 
von Neumann relative entropy introduced in |, | as a measure of entanglement. Starting from 
BD states, we perform local quantum operations and classical communications (LQCC) and find 
L-S decomposition for a generic two qubit system. We prove that for some special LQCC the 
obtained decomposition is optimal. It is also shown that for these cases average concurrence of the 
decomposition is equal to concurrence. 

The paper is organized as follows. In section 2 we review BD states and present a perspective 
of their geometry. L-S decomposition of these states is obtained in section 3 via a geometric 
approach. We prove that this decomposition is optimal. Relation between L-S decomposition and 
relative entropy is discussed in section 4. It is shown that BSA also minimize von Neumann relative 
entropy. Effect of LQCC on L-S decomposition is studied in section 5. The paper is ended with a 
brief conclusion. 

2 Bell Decomposable States 

In this section we briefly review Bell decomposable (BD) states and some of their properties. A 
BD state is defined by 

4 4 

p = X>lV>i)(^l, 0< Pi <l, 5>* = 1, (2-2) 

8=1 8=1 
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1 

0+> = _(|tT) + |U)), 


(2-3) 




r> = ^(ITT)-|U)) ; 


(2-4) 


l^3> = 


tf + > = ^(IU> + UT», 


(2-5) 


l^ 4 ) = 


iT> = ^(IU>-Ut»- 


(2-6) 



In terms of Pauli's matrices, p can be written as 



1 3 

p = -{I ®I + ^2t i a i ®a i ), (2-7) 
i=i 



where 



*i = P1-P2+P3-PA, 

*2 = "Pi + P2 + P3 - P4, ( 2 " 8 ) 
*3 = Pi + P2 ~ P3 ~ P4,- 

From positivity of p we get 

1 + h - t 2 + t 3 > 0, 
1 - ti + h + t 3 > 0, 

(2-9) 

1 + ti + t 2 - t 3 > 0, 

1 - ti - 1 2 - 1 3 > 0. 

These equations form a tetrahedral with its vertices located at (1,-1,1), (—1,1,1), (1,1,-1), 
(—1, —1, —1) [|10|| . In fact these vertices are Bell states given in Eqs. (|2-3| ) to (|2-6j ), respectively. 



According to the Peres and Horodecki's condition for separability [11, 12], a 2-qubit state is 



separable if and only if its partial transpose is positive. This implies that p given in Eq. (2-7) is 
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separable if and only if tjS satisfy Eqs. fl2-9|) and 



1 + h + t 2 + t 3 > 







1 - h - t 2 + t 3 > 







(2-10) 



1 + t x - t 2 - i 3 > 







i-t l + t 2 -h> 



o. 



Inequalities ( |2-9| ) and ( 2-10| ) form an octahedral with its vertices located at Of 1 = (±1,0,0), 



O2 = (0, ±1,0) and = (0,0, ±1). Hence, tetrahedral of Eqs. (2-9) is divided into five regions. 



Central regions, defined by octahedral, are separable states. There are also four smaller equivalent 
tetrahedral corresponding to entangled states. Each tetrahedral takes one Bell state as one of its 
vertices. Three other vertices of each tetrahedral form a triangle which is its common face with the 
octahedral (See Fig. 1). 

3 Lewenstein-Sanpera decomposition 

According to Lewenstein-Sanpera decomposition Q, any 2-qubit density matrix p can be written 



where p sep is a separable density matrix and \tp) is a pure entangled state. The L-S decomposition of 
a given density matrix p is not unique and, in general, there is a continuum set of L-S decomposition 
to choose from. The optimal decomposition is unique in which A is maximal, and 



as 



P = ^Psep + (1 — A) (V| 



A G [0, 1] 



(3-11) 





(3-12) 



All other decomposition of the form p = \p sep ± (1 — A) ijjj yj) , with A G [0, 1] such that p 7^ 
necessarily implies that A < A^ op ^ Q. 
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In the following, we will refer to Eq. ( |3-11 ) as the optimal decomposition of p. The separable 



part p sep is called the best separable approximation (BSA) of p, and A is its separability. 

Here in this section we obtain L-S decomposition for Bell decomposable states via a geometrical 
approach. Our results are in agreement with those reported by [Q, ||. In addition we present an 
explicit form for p sep and show that, pure entangled state is Bell state which p belongs to its 
entangled tetrahedral. For simplicity, we show in Fig. 2 entangled tetrahedral corresponding to 



singlet state (2-6). 



Suppose p is an entangled state parameterized as t = (ti,t2,t 3 ). We connect vertex p, which 
denotes singlet state, to point ~t and extend it to cut separable surface O^O^O^ at 7 correspond- 
ing to separable state p s , where this line can defined by Eqs. (l + ^X^i — *i) — (1 + ti ) (a?2 — £2) = 
and (1 + t 3 )(x2 — £2) — (1 + £2X^3 — £3) = 0. It can be easily seen that this line cuts plane O^O^O^, 
defined by x\ + X2 + £3 + 1 = 0, at point t = (t'i,t' 2 ,t' 3 ) 

t , = -i + h-t 2 -t 3 

1 3 + h + t2 + t 3 ' 

2 3 + ii + t 2 + *3 ' 

t , = -1-ti -t 2 + t 3 

3 3 + h + t2 + t 3 ' 

Using Eq. ( |2-8| ) it is straightforward to obtain coordinates of p s in terms of parameters pi as 

Pi= 2(1 -p 4 ) f ° r i = 1 ' 2 ' 3 and & = \ ( 3 " 14 ) 
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Now, using the Eqs. (l^T)) , (|^D we can write explicit form for separable state p s 



( 



1 



2(3 + ti +t 2 + t 3 ) 



l + t 3 h-t 2 

2 + h + t 2 -l-t 3 

-l-t 3 2 + t 1 +t 2 

h-t 2 o o i + 1 3 



(3-15) 



By convexity we can write p as convex sum of p s and projector \ip ) {ip 



p = \p s + (l-\)\^-){^-\. 



(3-16) 



Using Eqs. ( [3-13 ) and ( 3-1 6| ) we obtain 



A 



\pt\ _ 3 + h + t 2 + t 3 

\pt'\ ~ 2 



1-C, 



(3-17) 



where \pt\ and \pt'\ are distances between points p, t and also p, t' respectively, and C is concurrence 

Obviously Eq. ( |3-17D implies that the entanglement contribution of singlet state in L-S decom- 
position of the BD states is the same as its concurrence. The concurrence of a mixed state is defined 
as the minimum of the average concurrence over all decompositions of the state in terms of pure 
states pL This means that for the L-S decomposition given in (|3-11|) we have C(p) < (1 — A)C(^). 



Eq. (3-17) shows that optimal decomposition of BD states saturate this inequality. 

Now, in order to show that A of Eq. (3-17) is maximal and thus the decomposition ( [3-16] ) with 
p s given in Eq. fl3-15| ) is optimal, first we show that p s can be written in terms of product states. In 
fact p s of Eq. (|3^1~5|) can be written as a convex sum of three states corresponding to three vertices 
Oi , O2 , O3 of octahedral, 

p s = \{p{ + \ 2 p 2 + \ 3 p 3 , (3-18) 
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where 



p 1 = t;(\x + ) (x+\ <8> \x-) + 8) (x+|), 

p2= U\y+)(y+\®\y-)(y-\ + \y-) (y-\®\y+) (y+\)> (3 ~ 19) 

P3 = U\ Z +) ( Z +\ ® + ® 

and |x±), |y±) and \z±) are eigenstates corresponding to eigenvalues ±1 of cr x , a y and a z , re- 
spectively. Considering the fact that for any point interior to equilateral triangles, the sum of 
the orthogonal distances of the point to three edges is equal to the height of triangle, for trian- 
gle O^O^O^ we have hi + h 2 + hs = \/3/2, where hi, h% and h% are orthogonal distances from 
corresponding edges (See Fig. 3). After straightforward calculation we get 

h - - ft *> 
V 2 



where t^, (i = 1,2,3) are coordinates of p s given by ( |3-13 ) 



Taking into account the fact that \ is proportional to hi and A = A 1 + A 2 + A 3 , we get 

A^ = |(l-tl +t 2 + t 3 ), 

X -= i ( i + tl _ t2 + t3)) (3-20) 

Ag:= l(i + tl + t 2 -t 3 ). 
The same is true for other Bell decomposable states belonging to other maximally entangled 

tetrahedral. 

Using above results we rewrite p given in Eq. ( |3-16 ) in terms of product states and pure 
entangled state 

6 

P = ^ ' A a \ e a , f a )(e a J a \ + \i>-)(^-\, (3-21) 

a=l 
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where \e a , f a ), a = 1, 2, 6 are product states defined by 

|ei,/i) = \x + ) <g> \x-) , |e 2 , /2> = <8> \x+) , 

|e 3 ,/s) = b+>®|l/-) , |e4,/4> = |l/-)®b+>. ^ 3 " 22 ^ 

les, /s) = |*+) ® k-> , |e6,/e) = \z+) ® |-Z-) , 
and A a , a = 1, 2, ...6 are given by 



Ai = A 2 = (3-23) 

A 3 = A 4 = (3-24) 
K 



A 5 = A 6 = '-j-. (3-25) 

Now with this notation, we are in position to prove that decomposition ( |3-16| ) with A given in (|3] 
17|J is optimal. To do this we show that all coefficients of product states appeared in (3-21) are max- 



imal. According to 0] maximizing all the pairs (A Q , A^) with respect to p a p = p— J2 a '^a,/3 Ay -Fa' 
and (PajPg) is a necessary and sufficient condition to subtract the maximal separable matrix 
from p, where for the sake of self-containty we quote the theorem 2 and the related lemmas of 
reference Q below. 

Theorem 1 

Given the set Ay of product vectors \e, f) £ Tt{p) , the matrix p* = J2a^-aPa is the best 
separable approximation to p iff a) all A a are maximal with respect to p a = p — J2 a i^ a A Q /P a ' and 
the projector P a ; b) all pairs (A a , A^) are maximal with respect to p a p = p — J2a'^a,/3 ^-a'Pa', and 
the projectors (P a ,Pp). 

Lemma 2 A is maximal with respect to p and P = \ip) (ip \ iff a) if \yb) 7Z(p) then A = 0, and 
b) if \ip) e K(p) then A = p~ x > 0. 
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Lemma 3 ^ A pair (Ai, A2) is maximal with respect to p and a pair of projectors (P\,P2) iff: a) 
if\ipi), \ip2) do not belong to TZ(p) then Ai = A2 = 0; b) if does not belong, while |V>2) G ~R<(p) 
then Ai = 0, A 2 = (ip 2 \ P~ X IV^) -1 / c) if \rb 2 ) G K(p) and (ipx\ p' 1 \ip 2 ) = then A { = 

{tpi\ p' 1 \ipi)~ , i = 1,2; d) finally, if \ipi) , \ip 2 ) G ^(p) and (^ x | |^ 2 ) / i/ien 

Ai= ((H p- 1 M - \ (H p- 1 M \)/D, 

(3-26) 

a 2 = ((^ip- 1 !^)-! (Hp- 1 ^) \)/d, 

where D = (V>i| p" 1 (^| P' 1 \ih) ~ I (^i|p-i I ^2) ? ■ 

First, we show that A a s are maximal with respect to p a and P a . 

Matrices p a = p-J2a'^ a = A a P a +(l-\) \ip~) (ip~ \ with P a = \e a , f a ) {e a , f a \ for (i = 1,2,. ..,6) 
have two zero eigenvalues and two non zero eigenvalues. In Bell basis its kernel and range are 
separated. After restriction to its range, it is straightforward to evaluate p^ 1 and we find that 

((H,fi\ Pi 1 \(%,fi) = 1/Aj. 

In order to prove that the pair (A a ,A^) are maximal with respect to p a p and the pair of 
projectors (P a ,Pp), we proceed as follows: 

a) Matrices p%,i+x = AjPj + Aj+iPj+i + (1 — A) {tp~\ for (i = 1, 3, 5) have a two dimensional 
range . In Bell basis its range and kernel are separated and one can obtain (e^, /j| p^j+i \ e i: fi) = 
(Aj+i + (1 - X))/Ti, (e i+1 , f i+1 \ p~rj +1 \e i+1 J i+1 ) = (A { + (1 - A))/^ and (e,, ft\ p~j +1 \e i+1 , = 
(1 — A))/(2Fj), where Tj = AjAj + i + |(1 — A). Using the above results together with Eqs. ( p-26| ) we 
obtain the maximality of pair (A,, Aj+i) with respect to Pi^+i and the pair of projectors (Pi, Pj+i) 
for i = 1,3 and 5. 

b) For other possibility of a and (3, matrices p a p = A a P a + ApPp + (1 — A) (Y>~ | have rank 
3 . Using the Bell basis we can evaluate p^ 1 and we find that {e a , f a \ p _1 \e@, fp) = for a ^ /?, 
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(ea, fa\ p 1 \e a , f a ) = 1/A a . This completes the proof that A Q of Eq. ( 3-20 ) are maximal and 
decomposition ( |3-1£| ) is optimal. 



Also it is worth to note that the decomposition ( 3-16| ) satisfies conditions for BSA of Ref. 



According to theorem 1 of Ref. Q, decomposition given in ( |3- 16 ) is the optimal decomposition if 
and only if: rank( / oJ s ) = 3, i.e. 3u\ p^ B \4>) = 0, and either 

(i) 3 Q>0 (!</>) (<A|) Ts IV) = "« IV), or 

(3-27) 

(ii) rank(p s ) = 3, i.e. 3^ Ps \~4>) = 0, and 3 a) „> (y\4)Q>\ + (\<t>) (<t>\f B ) IV) = -*W) ■ 

It is now straightforward to see that p T s B has three non vanishing eigenvalues, that is, its rank is 
3. Its one dimensional kernel is along the Bell state |Vi) given in Eq. ( |2-3j ) . Actually the density 
matrices corresponding to the interior of tetrahedral satisfy condition (i) while those at its boundary 
satisfy condition (ii), respectively. 

4 Relative entropy of entanglement and L-S decomposition 

Vedral et al. in |9| introduced a class of distance measures suitable for entanglement measures. 
According to their methods, entanglement measure for a given state p is defined as 

E(p)=xmnD(p\\a), (4-28) 

where D is any measure of distance (not necessarily a metric) between two density matrix p and a, 
and T> is the set of all separable states. They have also shown that von Neumann relative entropy 
defined by 

S(p || a) = tr{p In-}, (4-29) 
a 

satisfies three conditions that a good measure of entanglement must satisfy ||. Here, we would like 
to emphasis that p s given in Eq. ( |3-15| ) minimizes von Neumann relative entropy given in ( 4-29| ). 
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Authors in ||] have shown that for BD states given in Eq. (2-2), separable state a that minimize 
relative entropy is 

Pi = on ~ s for i = 1,2,3 and p' i = \. (4-30) 
2(1 -pi) 2 



It is worth to note that the above equation is the same as Eq. ( 4-30| ), that is, separable state 



optimizing L-S decomposition minimizes von Neumann relative entropy, too. 



5 L-S decomposition under LQCC 



In this section we study the behavior of L-S decomposition under local quantum operations and 



classical communications (LQCC). A general LQCC is defined by [14 

, _ (A ® B)p(A ® B)i 
9 ~ tr{{A®B)p{A®B)\y 

where operators A and B can be written as 



(5-31) 



A®B = U A F> a ' m ®U B f u > b > n , 



(5-32) 



where U A and Ub are unitary operators acting on subsystems A and B, respectively and the 
filtration / defined by 

fn,a,m = fJi (i 2 + am.a), 



(5-33) 



As it is shown in Refs. 14, 15|, the concurrence of the state p transforms under LQCC of the form 



given in Eq. ( 5-31 ) as 



r< >\- M 2 ^ 2 (i-« 2 )(i-^ 2 ) r( , 

[P> tr({A®B)p(A®B) J t) {P> ' 
Performing LQCC on L-S decomposition of BD states we get 



(5-34) 
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tr((A <g> B)p{A ® B)i) 
with p' s and |^') defined as 



{ A®BWA®B)1 _ x , p , H1 _ xl)Wwl (5 . 35) 



A = %*m$*2L , (5-36) 



(^0B)|V>-) 



W) = iA*wn (5 _ 37) 



respectively, and A' is 



x , = tr { (A®B) Ps (A®B)l) _ (5 _ 38) 



tr((A® B)p(A(g> B)^ 

Using Eq. ( |5-3g| ), we get for the weight of entangled part in the decomposition (|5-35| 



n \'\ {^{AA^BB^-) 



Now we can easily evaluate the average concurrence of p' in the L-S decomposition given in ( 5-35| ) 



where, by Comparing the above equation with Eq. ( |5-34 ) we see that (1 — A)C(|V>}) (the average 



concurrence in the L-S decomposition) transforms like concurrence under LQCC. 

In order to prove that the decomposition ( 5-35| ) is the optimal one, we rewrite p s in terms of 



the product states given in Eq. ( |3-22 ) 



Ps = ^2 ^a\e a , f a ) {e a , fa\ ■ (5-41) 



Now, performing LQCC action we get 

6 



Ps = T, A '<xKf'«)KJ'«\, M2) 



a=l 

where 



lC f a ) = (A0 g^.' /a) , t (P a ) = (e a , f a \ (Arf ® B5t) |e a , / a ) , (5-43) 
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with P a = \e a , f a ) (e a , f a \ and 

tr{(A®B)p{A®B)i) a v ' 

First we show that A^s are maximal with respect to p' a and the projector P' a . 

As we see the matrices p a = A. a P a + (1 — A) \ip~) {ip~\ for (a = 1, 2, 6) transform as 

, (A®B)p a (A®By 



Pa 



t(p a ) = tr((A ® B)p a (A ® B)t), (5-45) 



under LQCC. Using the fact that LQCC transformations are invertible [|l5|, [H], [17|] , we can evaluate 
Pa as 

Pa 1 = *(Pa) (At ® S^-Va^A ® (5-46) 
Using the above equation and Eq. ( |5-43 ) we get 



( e ai fa\ Pa \ e ai fa) ~ ^rp \ ( e «» fa\ Pa \ e a, fa) — A a . (5 47) 



The Eq. ( 5-47| ) shows that A^,s are maximal with respect to p' a and the projector P' a . 



In order to prove that the pair (A' Q ,A^) are maximal with respect to {p' a , p'p) and (P a ,Pp), we 
proceed as follows: 

a) Matrices p^j+i = AjPj + Aj + iPj + i + (1 — A) \ip~) {ip~\ transform under LQCC as 

pWi = ^Mfcil^iM)! , t(p . ) = tr{{A ® B Ki+l(i4 ® £)t) 5 (5-48) 
HPi,i+lJ 

Using the above equation and invertibility of LQCC we arrive at the following results 

t(p) ( A i+i + |(l-A) 



( e j> fi | Pi,i+1 | e i) /i) 



t(Pi 



(e> f \r W f N-^) ( A< + g (1 - A) ) (5 49) 

\ e i+l' Pi,i+1 l e i+l> — i(P) p. ' 



L-S decomposition for BD states 



16 



{ e ii fi | Pi,i+1 | e i+l> fi+l) 



t(p) (1 - A) 



y/t{Pi)t(P i+1 ) 2Ti 

where Tj = AjAj + i + |(1 — A)(AjAj+i). Now using Eqs. (|3-26|) we get 

( e i) /il Pi,i+1 l e i) /i) — I ( e i-> fi \ Pi,i+1 | e i+li fi+l) I _ »/ 



D 



H+l 



(t(p i+1 ) - ^t(Pi)t(p i+1 : 

X " A) TPS 

Hp) 



+ (1-A) 



(5-50) 



( e i+l) /i+ll Pi,i+1 K+l; fi+l) I ( e i> /il Pi,i+1 | e i+l' /i+l) 



A' 



+ (1 " AP T7 \ 

hp) 



(5-51) 



b) For other values of a and (3, matrices p Q)i g = A a P a + ApPp + (1 — A) \ip ) (if) | transform 
under LQCC as 



Pa,f3 = —„ 7\ ' HPa,/?) = tr (( A ® B)p at p{A (g> BY). 



t{Pa,p) 



So with same procedure we can evaluate the following expressions 



(5-52) 



( e a' fa | Pa, (3 | e a> fa) 



e 'pJf3 



Pa,P 



A' ' 



(5-53) 



<e' Q , /al Pai e /3' /s) = °' for a ^ P- 



Eqs. ( |5-50 ) and ( 5-51 ) show that the pair (A' a ,A'a) are maximal with respect to p a ^ and 
(P^,PL) provided that t(Pi) = t(Pi + i). This restricts LQCC to special case of A = B. Under these 
conditions the decomposition given in Eq. ( |5-35 ) is optimal. Actually considering the fact that the 
inequality C(p') < (1 — X')C(ip') is saturated by decomposition given in (|5-35 ) we conjecture that 
this decomposition is optimal for a general LQCC, that is, p' s is the BSA for p' . This implies that, 
in general, the product states given in Eq. ( |5-43| ) is not the good product ensemble for the best 
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separable part. To prove that the decomposition ( 5-35| ) is the optimal one for the whole class of 
LQCC is still an open problem which is under investigation. 



6 Conclusion 

We have derived Lewenstein-Sanpera decomposition for Bell decomposable states from an entirely 
different approach. We show that for these systems, the weights of the pure entangled part in the 
decomposition is equal to the concurrence of the states. Optimality of the presented decomposition 
have been proved by using the theorems given in pi] . It is also shown that the optimized separable 
part of L-S decomposition minimizes the von Neumann relative entropy. We have also obtained 
Lewenstein-Sanpera decomposition for a large class of states obtained from BD states via some 
LQCC action. It is shown that for these states the average concurrence of the decomposition is 
equal to their concurrence. 
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Figure 1: All BD states are denned as points interior to tetrahedral. Vertices P\, P2, P3 and 
P4 denote projectors corresponding to Bell states Eqs. (2-3) to ( p-6|) , respectively. Octahedral 
corresponds to separable states. 

Figure 2: Entangled tetrahedral corresponding to singlet state. Line pc denotes entangled 
Werner states. Points t and t' correspond to a generic BD state p and associated BSA p s . Vertex 
p denote singlet state and other vertices are defined in Eq. ( 3-19| ). 



Figure 3: p s can be written as a convex combination of separable states p 1 , p 2 and p 3 with 
weights proportional to hi, hi and /13, respectively. 




Figure 1: 




Figure 2: 




Figure 3: 



